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Continuous time random walk (CTRW) subdiffusion along with the associated fractional Fokker- 
Planck equation (FFPE) is traditionally based on the premise of random clock with divergent mean 
period. This work considers an alternative CTRW and FFPE description which is featured by finite 
mean residence times (MRTs) in any spatial domain of finite size. Transient subdiffusive transport 
can occur on a very large time scale r c which can greatly exceed mean residence time in any trap, 
t c 3> (t), and even not being related to it. Asymptotically, on a macroscale transport becomes 
normal for t ^> r c . However, mesoscopic transport is anomalous. Differently from viscoelastic 
subdiffusion no long-range anti-correlations among position increments are required. Moreover, our 
study makes it obvious that the transient subdiffusion and transport are faster than one expects 
from their normal asymptotic limit on a macroscale. This observation has profound implications 
for anomalous mesoscopic transport processes in biological cells because of macroscopic viscosity of 
cytoplasm is finite. 

PACS numbers: 02.50.Ey, 05.40.Fb, 05.60.Cd 



The subject of anomalous diffusion and transport be- 
came increasingly popular in the last two decades (see 
Fig. 1 in [lj]) with many application fields in physics 
and beyond and many theoretical approaches coexist- 
ing and developed in parallel. The works by Montroll 
and Weiss 0, Scher and Montroll and Shlesinger Q 
laid the grounds for the continuous time random walk 
(CTRW) approach to subdiffusion based on the premise 
of a random distribution of the residence times spent in 
traps with non-existing mean time period between two 
successive spatial steps. Assuming that the residence 
times in traps are independently distributed and no cor- 
relations between spatial steps exist (semi-Markovian as- 
sumption) , the whole process is completely characterized 
by a residence time distribution ^i{ T ) m the i-th trap 
and the probability distribution p(k\i) to jump from i 
to another trap k. The highly dispersive subdiffusive 
transport with the position variance growing sublinearly, 
(5x 2 (t)) oc n a t a , < a < 1, where is subdiffusion co- 
efficient, emerges when the residence time distributions 
behave asymptotically as 



t/»(r) oc r a ~ 



(1) 



This approach acquired enormous popularity and has 
been followed by numerous workers in the field 043- 
However, there are no real processes in nature which 
would show this characteristic feature, i.e. infinite mean 
residence time (MRT) in a finite spatial domain. Of 
course, in reality it is always finite. Nevertheless, it can 
largely exceed a characteristic time t l ~ (L 2 /K a ) 1 / a re- 
lated to a finite linear size L of the medium, e.g. the 
width of a thin amorphous semiconductor film exhibit- 
ing akin transport properties. Then, the approxima- 
tion of infinite mean residence time can be well justi- 
fied. Asymptotically transport will be normal. How- 



ever, this macroscopic normal transport regime is irrele- 
vant on micro and nanoscales in the corresponding amor- 
phous media. The main advance of this work is to in- 
troduce an alternative semi-Markovian CTRW approach 
to subdiffusive transport featured but by finite MRT as 
a main premise. Moreover, we derive a modified frac- 
tional Fokker-Planck equation (MFFPE) in the continu- 
ous space limit, where MRT in a point-like trap becomes 
infinitely small. Within this alternative CTRW approach 
MRT (r) in any finite spatial domain not only remains 
finite, but it shrinks to zero with vanishing size of the 
domain. Nevertheless, the corresponding residence time 
distribution does scale as one in Q on a long interme- 
diate time scale. Since in all real experiments the mean 
residence times are always finite one can conjecture that 
some experimental observations of the law ([1]) along with 
subdiffusion, e.g. in [8], might be more consistent with 
our alternative description, rather than with the popular 
one. 

The key modification in this alternative is to consider 
the residence time distribution on sites characterized by 
the Laplace-transformed survival probability 



i 



s + L, 7 =o r " 



(2) 



where a.Q = 1 and < ay < 1 for j — 1, 2, N. By not- 
ing that ipi(s) — 1 — s$i(s), the corresponding residence 
time distribution ipi (r) can be represented as a sum over 
N + 1 escape channels, ipi{r) = J2i=o ^j^ 1 "); with the 



Laplace-transformed ipji( T ) reading 



(3) 



This justifies (see also below) to name the quantities r .j 



the fractional rates. One of them, r*i ; 



3) i is just the 



2 



normal rate corresponding to Qo = 1- The presence of 
such a normal escape route from any trap is a crucial 
point yielding a finite MRT, (n) = $i(0) = 1/7%; in 
the trap i. This is because the particle can go cither 
via any slow route, or through the fast one which de- 
fines the MRT. The MRT is thus finite in our model, but 
its variance is not bounded, (Srf) = oo, reflecting the 
presence of anomalously slow pathways. Similar hap- 
pens e.g. if the particle can escape out of the anomalous 
CTRW subdiffusion domain being injected on the nor- 
mal radiation boundary which it can either leave with 
normal rate immediately or penetrate further inside the 
anomalous domain ||. Even more strikingly, by doing 
the continuous space limit below we shall consider the 
limit (Tj) — > 0. Here lies the principal difference of our 
subdiffusive CTRW model from a popular one dominat- 
ing currently over this research domain. The simplest 
representative of this class of models is (N = 1) 



(n) 



l + [T c s) 1 - a + s(n) ' 



(4) 



where t c is a cutoff time which is equal for all traps and 
defines the maximal time range of subdiffusion. For a = 
1/2 this expression can be easily inverted to the time 
domain and expressed in terms of a combination of the 
error function and power law functions. Generally, it 
is not difficult to show that for (n) <C r <C T c , ipi(r) 
exhibits precisely the same scaling behavior as one in 
Eq. (fT]). Asymptotically for r 3> t c it is but a different 
power law, ipi(r) oc t 3 ~ q ensuring that the MRT is finite. 
Apart from this, initially for t <C (rj), ip(r) oc t~ 1+ ". 
All in all, our RTD exhibits three different interchanging 
power-laws even in the simplest case of N = 1, providing 
a rather rich model. 

Next, let us consider the generalized master equation 
(GME) for the CTRW on the lattice x t = iAx (i e Z) 
with a discretization step Ax and nearest neighbor jumps 
only, characterized by the jump probabilities qf — p(i + 
and q~ — p(i — l\i) obeying qf + q[ = 1. The 
corresponding GME reads [10( 



{K+^t - i')^-i(t') + Kr +1 {t - t')P l+l {t') 



- [Ktit-t^ + Krit-t'mit')}^', 



(5) 



where the Laplace-transform of the kernel Kf(t) is re- 
lated to the Laplace-transform of the residence time dis- 
tribution (RTD) by Kf (s) = qfip,(s)/^(s), which for 
the considered model yields 



N 



Kt{s) = qf ( /•:, • X''" 



(6) 



where qf and q i are the forward and backward jump 
probabilities, accordingly, q^~ + q^ = 1. One can define 



also forward an backward fractional rates so that r 
lt r aj,i and r a]ti = r+ ^ + r"^, and rewrite the 
([5]) in the form of a fractional master equation 



N 



+r-.. +1 P i+1 {t) - [r^i + r-jPiit) 1 , (7) 



where 



t b]f{t)~ ' 



'dif J^L , (8) 



r(i- 7 )^y t0 (t-f)i" 

with < 7 < 1 defines the fractional Riemann-Liouville 
derivative acting on some function f(t) Q. 

Furthermore, let us consider an additional potential 
field V(x) which influences the jumping rates as follows 

r%,i = {K a J{Axf)e*. V [-p{V i±l/2 -V{)] (9) 

where n aj has the meaning of a (generally) anomalous 
diffusion coefficient, and (3 = 1/ (ksT) is inverse temper- 
ature. The form ^ ensures that the Boltzmann equilib- 
rium condition, r+ j . i _ 1 /r~. l = exp[(3(Vi~i - Vi)], is sat- 
isfied for any V(x) and our description is thermodynami- 
cally consistent. Using a finite difference operator A/ Ax, 
AP(x,t)/Ax := [P{x + Ax/2,t) - P(x - Ax/2,t)]/Ax, 
one can rewrite the GME ([5]) for the considered model as 



P(Xi,t) = 





N 




+ K "3 




J'=l 


-( 


' ' e -PV{xi) 


Ax \ 





't 



where we have explicitely singled out the normal transi- 
tion channel which is present between any two localized 
states and acting in parallel to all other ultraslow chan- 
nels. 

In the continuous space limit Ax — > with constant 
K aj Eq. (fTU)) yields the following generalization of the 
fractional Fokker-Planck equation (FFPE) of Refs. 0, 
11 13] 




w, 'P(i,t) . (11) 



dP(x,t) 

Ft 



The generalization consists primarily in the presence also 
of a normal diffusion term along with the anomalous dif- 
fusion contributions. Notice that it was not only obtained 
from the CTRW possessing finite MRT in each trap but 
even in the limit (r,-) oc (Ax) 2 /ki — > 0, so that the MRT 



in any spatial domain of a finite size is obviously finite. 
For N = 1, and in the limit K\ — > it reduces to the 
original FFPE of works llj, |l2j. In the same limit, but 



for N = 2 it reduces to the modified FFPE of Ref. [14 1. 
Generally, it belongs to the class of distributed order frac- 
tional kinetic equations 1J, [l5( . 

In the following, we consider the simplest case of 
N = 1, with finite k,\ and n a — KiT^" with a very 
large r c . The key point is that for t <C t c the normal 
diffusion term is negligible and the dynamics is subd- 
iffusive, whereas the mean residence time in any finite 
spatial domain is finite. Here lies a profound, striking 
difference with the earlier approach. For t 3> t c , the 
normal diffusion contribution dominates the transport, 
which becomes asymptotically normal. However, the bi- 
ased subdiffusion changes asymptotically into superdiffu- 
sion (see below) because the variance of residence times 
is infinite. 

Our focus is now on an intermediate subdiffusive 
regime. Considering the mean displacement under the 
constant force F, the solution of (fTT|) yields 



(6x(t)) = Ft/m + Fe/i^Til + a)] 

t 

77 



FTc 

m 



i 



r(l + a) 



(12) 



where the normal, 771, and fractional, r] a , frictional co- 
efficients are related to the corresponding diffusion co- 
efficients by standard, 771 = hsT/ni, and generalized, 
T\a = ksT/Ka, Einstein relations, respectively. Fig. [T] 
displays this solution for three different values of a along 
with one and the same normal diffusion asymptotics for 
t t c and three different initial t <C t c asymptotics. One 
can conclude that the motion remains subdiffusive for 
most time and can be characterized by a time-dependent 
power-exponent a e s(t) which gradually changes from a 
for t <C t c to a = 1 for t 3> r c . Even more strikingly, 
in the subdiffusive regime the transport is as a mat- 
ter of fact faster than its normal diffusion asymptotics. 
Indeed, already the structure of the generalized master 
equation ([5]) says that in the Markovian approximation, 
Pj(t') w Pi(t), only the normal rate term contributes 
asymptotically, t — > 00. This is because the terms in the 
integral of the memory kernel, J*Kf(t')dt', which corre- 
spond to the anomalous rate contributions, tend to zero 
with increasing the upper limit t, as l/t 1 ~ aj , being al- 
ways positive. Initially these terms, however, profoundly 
contribute making the overall transport faster than its 
normal asymptotics. Notice that this interpretation is 
opposite to one suggested in Ref. [15] for diffusion which 
we consider shortly: the transport obviously slows down 
in the course of time and not accelerates. 

For the displacement variance (6x 2 (t)) = (x 2 (t)) — 
{5x(t)) 2 , which characterizes diffusion, after some alge- 
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FIG. 1: (Color online) Anomalous transport which asymptot- 
ically becomes normal. 



bra we obtain 



(Sx 2 (t)} 



2K 1 t + 2K a t a /r(l + a) 

m 
2 



(13) 




l + Q 



r(i + 2a) r 2 (i 



2a 



Taking into account the generalized Einstein relation, it 
is seen immediately that a linear fluctuation-response re- 
lation, (Sx(t)) = [F/(2k B T)](Sx 2 (t)) F=0 is satisfied, be- 
tween the "fluctuation" (dx 2 (t))p = o in the absence of 
driving force F, and the "response" (5x(t)) to it. In the 
subdiffusive transport regime restricted by t <C t c , 



(5x(t)) w Ft a /[r) a T(l + a)}, 
(5x 2 {t)) w 2K a t a /T(l + a), 



(14) 



i.e. the variance (6x 2 (t)) behaves subdiffusively also in 
the presence of bias F, like the mean displacement in 
Fig. [TJ This is in a striking contrast with the standard 
CTRW and FFPE subdiffusion based on the premise 
of infinite mean residence times. In the latter case, 
{Sx 2 (t))F^o oc Ft 2a , for a sufficiently large F, or asymp- 
totically, similar to the last term in Eq. (|13p . Therefore, 
our theory is closer to what one intuitively expects for 
normal transport, sharing this normal feature also with 
the alternative Generalized Langevin Equation (GLE) 
approach to subdiffusion which we shall discuss shortly. 

Within the presented modified theory the large time 
asymptotic behavior of the variance in the presence of 
bias is also rather unusual. Obviously for F ^ 0, and 
t ^> r c , the main contribution in Eq. (|13p is superdiffu- 
sive, (5x 2 (t))p 9 Lo oc Ft 1+a . The explanation of this is the 
following. Even if the random clock governing subdiffu- 
sion in our case does possess a mean value, its variance 
is but infinite and in accordance with Ref. [lj (see Ap- 
pendix therein) determines the displacement variance in 
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the case of biased diffusion. Unbiased diffusion remains 
but asymptotically normal. Such a combination of nor- 
mal transport and superdiffusion was found recently in 
numerical simulations of normal Langevin dynamics in 
near critically biased washboard potentials in the pres- 
ence of frozen disorder for some parameters (l6j . In this 
respect it is worth mentioning that hopping diffusion in 
disordered potentials is equivalent within the mean-field 
(effective medium) approximation to a CTRW 5]. This 
provides one of physical justifications for the CTRW phe- 
nomenological approach. 

Let us compare now the considered modification of 
CTRW subdiffusional transport phenomenology with 
even more normal route of the GLE approach to sub- 
diffusion 17-2(1, where all the moments of the resi- 
dence time distribution in any finite-size spatial domain 
are finite 19]. Such a GLE approach is intrinsically 
based on the phenomenon of viscoelasticity [2l], [23} with 
long-range correlations in the Brownian particle displace- 
ments. Here, the physical origin of subdiffusion is en- 
tirely different from one in the semi-Markovian CTRW 
[T9I [20j . Stochastic dynamics is described by the over- 
damped gle [rzlGll 



rj(t - t')x(t')dt' 



dV(x,t) 
dx 



(15) 



where thermal force is (i) Gaussian, (ii) unbiased on 
average, and (iii) obeys the fluctuation-dissipation rela- 
tion, 



mat')) = kBTr ) (\t-t'\), 



(16) 



at the environmental temperature T. Macroscopic 
power-law models of viscoelasticity were introduced in 
1936 by A. Gemant 2l|, l22j as a generalization of the 
Maxwell model [23I ] with exponentially decaying vis- 
coelastic memory kernel. One of the Gemant 's macro- 
scopic models having relation to our present work corre- 
sponds to the memory friction 



fj(s) = 



1 + (sT C y- 



(17) 



in the Laplace space. For t <C t c , rj(t) can be approxi- 
mated by r)(t) ~ r) a /\T(l — a)t a ], with n a = rj\T^~ l and 
its asymptotic decay for t ^> r c is another power law, 
n(t) oc t a ~ 2 ensuring that 77(0) = rji is finite. This GLE 
model with the memory kernel (|17[) for f(x) = F = const 
yields again exactly the same solution (fT2"j) for the mean 
displacement and a modified (with F — > 0, the first 



line only remains) solution (fTB")) for the variance [2J|. It 
must be emphasized that GFFPE (fTTj) with N = 1 does 
not present a master equation counterpart of the GLE 
(fl~5|) with the memory kernel (fl7]) . as they have physi- 
cally (and mathematically too) totally different origins. 
The both approaches share, however, a prominent fea- 
ture that the mean residence time in any finite spatial 
domain is finite. 



This ensures, for example, that the asymptotic re- 
sponse of the corresponding anomalous dynamics to time- 
periodic perturbations is finite. However, these asymp- 
totical responses in both models can also be different in 
some details, reflecting different physics. For example, it 
has been recently shown that the traditional CTRW sub- 
diffusion based on infinite MRTs displays asymptotically 
dying response to periodic fields [25| which is at odds 
with most experimental data on anomalous dielectric re- 
sponse, and nevertheless agrees with some [26]. Nor- 
mally, the anomalous response exhibits aging behavior. 
However, the response normally decays not to zero, but 
to some finite value. The response to periodic fields can 
serve to distinguish among different underlying physical 
mechanisms. A generalization of the presented here mod- 
ified CTRW and FFPE model to time-dependent fields 
is not a trivial matter at all and it is left for a separate 
study. 

Furthermore, for subdiffusion and transport restricted 
by the time scale t < t c and spatial scale x < L ~ 
(kqt") 1 / 2 it is not easy to distinguish between two dif- 
ferent models discussed here which display a very similar 
transient behavior. Calculation of the velocity autocor- 
relation function as e.g. in Ref. 27| can clearly prove 
that the viscoelastic mechanism is at work, or statisti- 
cal analysis of single trajectories [28| can help to make 
the distinction. The common feature is, however, that 
in both cases the transient transport is faster than the 
asymptotically normal one, and the transient subdiffu- 
sion is faster than asymptotically normal diffusion in the 
unbiased case. The results on mRNA molecules subd- 
iffusion in bacterial cells 29j seem consistent with this 
paradoxical conclusion [24 1. Moreover, the experimental 
results on subdiffusion of the colloidal gold particles hav- 
ing a typical size of proteins with radius R = 2.5 nm in 
cytoplasm of living cells [30j delivers a further clear con- 
firmation. Indeed, in Ref. [301 ] experimental results were 
parameterized as 



(6x 2 (t)) =2Dt d 



t 

td 



(18) 



Comparison with our model yields r c = [T(l + 
a)W r f] 1/(1_a) > Ki=D and n a = Dr^- a . The exper- 
imental values of a, td and r s can be found for various 
cells in [3(| . The corresponding asymptotic normal diffu- 
sion coefficient is D = rn/(4 rp), where tq = 0.21 /j,m, is 
the microscope aperture size [30|. From this we can ob- 
tain «i, K a and t c , see in the Table I. For gold beads of 
this size in water Ki jWa t C r = 90/im 2 /sec. Therefore, from 
the experimental values of K\ in cytoplasm of selected eu- 
karyotic cells (possessing cytoskeleton) one can conclude 
that cytoplasm is about Ki )W ater/ti ~ 750 times more 
viscous than water and about two times less viscous than 
glycerol. On the maximal time scale of experimental data 
in Ref. [3(j, which is about 1 sec, the particles subdiffuse 
much faster than one would expect from the normal dif- 
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TABLE I: Parameters derived from experimental data 



Cell 


a 


TD, 

ms 


Ts, 

ms 


t c , sec 


«i, 

fj,m 2 /sec 


/1m 2 /sec" 


HeLa 


0.51 


91 


0.271 


30.3 


0.121 


0.644 


THLE 


0.48 


89 


0.329 


16.7 


0.124 


0.536 


HepG2 


0.52 


88 


0.767 


37.6 


0.125 


0.715 


HeLa + 500 


0.66 


92 


0.954 


2875 


0.119 


0.931 


mM sucrose 















fusion in such a viscous environment. This paradoxical 
and overlooked fact should be taken into account in any 
discussion of the influence of transient subdiffusion on 
living processes as it appears that subdiffusion can bring 
only advantages, contrary to some common expectations 

m 

The author is confident that the modified CTRW and 
FFPE theory of subdiffusive transport considered in this 
work presents a viable alternative to the traditional one. 
ft shares more similarities with viscoelastic subdiffusion 
and should be taken seriously as an ample modeling 
framework with attracting physical features and clear 
experimental relevance. 
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